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PRELIMINARY  REPORT  OF  NUMERICAL  SIMULATIONS  OF  INTERMEDIATE  WAVE¬ 
LENGTH  ExfiL  GRADIENT  DRIFT  INSTABILITY  IN  IONOSPHERIC  PLASMA  CLOUDS 

INTRODUCTION 

Barium  cloud  injection  experiments  [Rosenberg,  1971;  Davis  et  al . , 
1974]  have  provided  much  data  concerning  the  dynamical  evolution  of 
plasma  clouds  in  the  F  region  ionosphere.  The  characteristic  steepen¬ 
ing,  elongation,  and  strlation  of  barium  plasma  clouds  have  been  ex¬ 
plained  by  applying  the  ExB  gradient-drift  instability  [Simon.  1963] 
to  plasma  cloud  geometries  [Haerendel  et  al.,  1967;  Linson  and  Workman, 
1970;  Volk  and  Haerendel,  1971;  Perkins  et  al.,  1973].  Numerical 
solution  of  the  fundamental  fluid  equations  modelling  the  ExB  insta¬ 
bility  in  plasma  clouds  have  reproduced  not  only  the  large  scale  gross 
observational  features  of  plasma  cloud  evolution  [Zabusky  et  al . ,  1973; 
Lloyd  and  Haerendel,  1973;  Goldman  et  al . ,  1974;  Doles  et  al . ,  1976; 
Ossakow  et  al.,  1975,  1977]  but  also  the  spatial  power  spectra  [ Scanna- 
pieco  et  al. ,  1976]  and  minimum  scale  sizes  of  plasma  cloud  striations 
[McDonald  et  al.,  1978].  However,  the  goal  of  these  numerical  simula¬ 
tions  was  to  model  the  dynamics  of  the  entire  plasma  cloud  and  wave¬ 
lengths  less  than  about  a  kilometer  were  not  accurately  resolved.  But 
barium  cloud  experiments  sponsored  by  the  Defense  Nuclear  Agency  have 
shown  [Baker  and  Ulwick,  1978;  Kelley  et  al , ,  1979]  that  the  power 
spectra  of  the  plasma  cloud  striations  extends  to  wavelengths  on  the 
order  of  meters  to  tens  of  meters.  It  is  of  interest  to  numerically 
model  these  intermediate  wavelength  (100-1000  m)  irregularities  in  order 

to  compare  with  and  supplement  experimental  observations. 

Note:  Manuscript  aubmitted  October  fi,  1979. 
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In  this  preliminary  report  we  present  a  two-dimensional  numerical 


simulation  of  the  intermediate  wavelength  ExE  gradient-drift  insta¬ 
bility  applicable  to  local  unstable  regions  of  large  F  region  barium 
plasma  clouds.  The  results  of  these  simulations  will  be  shown  to  be 
consistent  with  experimental  and  theoretical  studies  of  striated  barium 
clouds  and,  in  many  respects,  similar  to  recent  numerical  simulations 
fKeskinen  et  al..  1979]  of  the  intermediate  wavelength  collisional 
Rayleigh-Taylor  instability  in  local  regions  of  the  bottomside  of  the 
equatorial  F  region  ionosphere. 

MODEL  EQUATIONS 

We  wish  to  model  two-dimensional  ExB  gradient-drift  processes  in 
local  unstable  regions  of  large  ionospheric  F  region  plasma  clouds. 

The  restriction  to  large  clouds  (large  Pedersen  conductivity  compared 
with  that  of  the  background  ionosphere)  permits  both  the  neglect  of 
the  cloud  interaction  with  the  background  ionosphere  (second  level) 
and  variations  of  cloud  density  and  potential  along  the  magnetic 
field  lines.  For  wavelengths  greater  than  the  ion-gyroradius  (approxi¬ 
mately  10  m  for  Ba+  in  the  twilight  F  region)  a  fluid  description  can 
be  used  which  equations  have  been  given  previously  [Perkins  et  al.. 
1973;  McDonald  et  al.,  1978;  Chaturvedi  and  Ossakow,  1979]. 

By  adopting  a  Cartesian  coordinate  system  (x,y,z)  with  the  geo¬ 
magnetic  field  B  in  the  z-directlon,  the  ambient  electric  field 
along  the  x-axis,  and  ignoring  to  lowest  order  electron  and  ion 
inertia,  we  can  write  after  transforming  to  the  cE  jfi/B2  frame 


where  n(x,y,t)  is  the  plasma  cloud  ion  density ,  cpi(xfy,t)  is  the 

electrostatic  potential  of  the  plasma  cloud,  K  1  *  K  *  +  1C  .  * 

r  e  en  ei  * 

K  ■  Q  /v  »  with  v  ,  v  .  the  electron  collision  frequencies  with 
ea  e  ea  en  ei 

neutrals  and  ions,  respectively,  and  Te  *  -  T  is  the  temperature. 

All  other  symbols  retain  their  conventional  meaning.  By  linearizing 
(1)  and  (2)  and  expressing  n  =  nQ  +  nj  with  nj,  ,  <*  exp 
[i(k  x  +  k  y)  +  y,  fc]»  k  *  B.  =  0,  k  L  >>1  we  find 


Vk  ■  '  <v«n  +  ',«l)(k2cs2>/ne”l 


where  L  1  =  (l/no) (3nQ/3y) ,  k2  *  k^2  +  k^2,  Cg2  ■  T/m We  note  that 
the  growth  rate  y^  maximizes  for  modes  perpendicular  to  the  density 
gradient  (k  -  k^)  while  the  modes  parallel  to  the  density  gradient 
(k  *  ky)  are  damped  by  cross-field  diffusion. 


Numerical  Simulations 


By  defining  n' (x,y)  -  n(x,y)/nQ(y) ,  cpj ,  (x,y)  *  q>j(x,y)/BL, 
x'  *  x/L,  y'  ■  y/L,  t'  ■  ct/L,  where  nQ(y)  will  be  defined  later, 
equations  (1)  and  (2)  can  be  written  in  dimensionless  form  as  follows 
(after  dropping  primes  for  clarity) 


3<p  3(j)  39 
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P2  n  3x  (4) 


where  gj  *  2T/eBK^L,  $2  *  E  /B,  03  “  T/eBL  are  dimensionless  constants. 

Equations  (3)  and  (4)  were  solved  numerically.  The  normalized 

plasma  cloud  density  n/nQ  in  (3)  was  advanced  in  time  using  a  variable 

time  step  flux-corrected  predictor-corrector  scheme  [Zalesak,  1979] 

which  is  basically  second  order  in  time  and  fourth  order  in  space. 

Equation  (4)  was  solved  for  the  cloud  potential  cpi(x,y)  at  each  time 

step  using  a  regridded  Chebychev  relaxation  method  [McDonald ,  1977]. 

The  results  to  be  presented  here  were  obtained  using  mesh  of  64  x  64 

grid  points  with  constant  spacing  Ax  *  Ay  ■  15  m.  Periodic  boundary 

conditions  were  imposed  on  n/nQ  ar.d  in  both  the  x  and  y  directions. 

The  parameters  used  in  these  simulations  are  typical  of  200  km  barium 

releases:  T  =  103  °K,  U g  ■  9  x  106  sec  l,  (Ba+)  “  36  sec  1 ,  B  ” 

0.5  G,  E  *  2mV/m,  v  .  ■  2  x  10 3  sec  3,  v  ■  150  sec  1  [Banks  and 
o  ei  en  v - 

Kockarts,  1973;  McDonald  et  al.,  1978]. 
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Three  different  computer  runs  were  made  using  different  values 

of  the  equilibrium  plasma  cloud  density  gradient  scale  length  L  *  3, 

6,  10  km  where  nQ(y)  =  Nq  (1  +  y/L),  Nq  constant,  is  a  steady  state 

(3/3t  ■  0)  solution  of  (1)  and  (2).  The  simulations  were  initialized 

with  a  two-dimensional  perturbation  of  the  form  [Chaturvedi  and 

Ossakow,  1979]  6n(x,y,t  =  0)/n.  =  A,  ,  sin  k._y  cosk.x  +  A2  osin  2k_y 

oi,i  y  x  »  j 

with  k  ■  k  =  2n/960  m  with  j  =  2  x  10  5  and  A2  o  “  2  x  10  6  (where 
x  y  *  « 

6n  s  n-n  ) .  We  will  now  present  the  important  nonlinear  results  of 
o 

these  simulations. 

Fig.  la  gives  an  isodensity  contour  plot  of  6n((x,y)/no  in  the 

x-y  plane  at  t  =  0  for  L  =  3  km.  The  initial  contours  describe  a 

sequence  of  enhancements  (<5n/n  >  0)  and  depletions  (6n/n  <  0) 

o  o 

arranged  in  a  checkerboard  fashion.  Fig.  lb  and  lc  show  the  evolution 

of  6n/n  at  t  =  1500  and  1800  sec,  respectively.  The  density  fluctua- 
o 

tion  contours  of  5n/n  at  t  *  2000  ser  are  displayed  in  Fig.  Id  where 

o 

some  elongation  in  the  y-direction  (ExB)  and  steepening  can  be  seen. 
Similar  contour  development  was  also  observed  for  the  other  two  plasma 
cloud  density  gradient  scale  lengths  L  =  6,  10  km  but  on  longer  time 
scales. 

Figs.  2a-b  show  representative  one-dimensional  power  spectra  both 
parallel  P(k  )  to  the  ExIJ  drift  and  perpendicular  P(k  )  in  the  non- 

y  x 

linear  late  time  regime  for  L  =  6  km.  These  power  spectra  are  defined 


as  follows 


P(k  )  * 

/  ^  i 

6n(k  ,k  )/n 

X 

J 

i y 

1  x  y  o 

p<y  - 

i 

6n(k  ,k  )/n 

1  x  y  o 

In  the  late  time  nonlinear  state  for  the  three  scale  lengths  L  =  3,  6, 
-n  -n 

10  km,  P(k  )  *  k  x  and  P(k  )  *  k  ^  with  n  =“  2-3  and  n  =*•  2-2.5 
x  x  y  y  x  y 

for  wavelengths  X  =*  70-960  m.  These  spectral  indices  are  consistent 
with  recent  in  situ  experimental  measurements  [Baker  and  Ulwick,  1978; 
Kelley  et  al.,  1979]  and  analytical  studies  [Chaturvedi  and  Ossakow, 
1979]  of  barium  cloud  striations. 

Finally,  Fig.  3  presents  the  time  history  of  modes  Aj  j  and  A2  o 

9  9 

for  L  *  10  km.  Initially  the  A2  o  along  the  I2xB  (density  gradi- 

9 

ent)  direction  damps  as  given  by  linear  theory.  As  the  linearly  un¬ 
stable  Aj  |  increases  it  nonlinearly  triggers  the  growth  of  A2  q.  At 

*  9 

late  times  these  modes  arrange  themselves  in  approximate  agreement  with 

the  nonlinear  amplitudes  predicted  by  Chaturvedi  and  Ossakow  [1979]. 

Similar  time  development  of  Aj  j  and  A2  y  are  observed  for  L  =  3,  6  km. 

*  » 

SUMMARY 

We  have  performed  preliminary  numerical  simulations  of  the  inter¬ 
mediate  wavelength  12*8  gradient-drift  instability  in  local  unstable 
regions  of  large  F  region  ionospheric  plasma  clouds  using  parameters 
typical  of  200  km  barium  releases.  For  barium  cloud  density  gradient 
scale  lengths  L  *3,  6,  10  km  we  find  that:  (1)  linear  unstable 
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inodes  saturate  by  nonlinear  excitation  of  linearly  damped  modes  along  the 
y  direction  parallel  to  ExB  drift  and  the  initial  cloud  density  gradient; 
(2)  in  the  nonlinear  well  developed  state  elongation  along  the  ExB  di¬ 
rection  is  seen;  (3)  the  one-dimensional  power  spectra  P(k  )  *  k  nx  and 
-n 

P(k  )  «  k  y  with  n  -  2-3  and  n  *  2-2.5  for  wavelengths  X  ■  70-960  m. 
y  y  x  y 

These  results  are  consistent  with  recent  barium  cloud  experimental  measure¬ 
ments  [Baker  and  Ulwick,  1978;  Kelley  et  al.,  1979];  support  the  nonlinear 
analytical  work  of  Chaturvedi  and  Ossakow  [1979];  and  shew  furthermore 
the  similarity  between  the  ExB^  gradient  drift  instability  in  ionospheric 
plasma  clouds  and  the  collisional  Rayleigh-Taylor  instability  in  the 
equatorial  F  region  ionosphere. 

Future  studies  are  planned  which  include  variation  of  the  para¬ 
meters  used  in  these  simulations,  addition  of  inertial  effects,  and 
inclusion  of  the  cloud  interaction  with  the  background  ionosphere 
(second  level). 
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gradient)  and  x  axis  horizontal  (perpendicular  toE  xB).  The  tiBk  marks  represent  grid  point 
locations.  ~° 


Fig.  2a  Fig.  2b 

-  One  dimensional  perpendicular  (a)  P(kx)  and  parallel  (b)  P(ky)  power  spectra 
md  ky,  respectively,  for  L  =  6  km  in  late  time  regime  (t  ■  3150  sec.).  Solid 
;  least  squares  fit  to  results  from  numerical  simulation  (dots)  and  Icq  *  2n/960 
fundamental  wave  number. 
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Fig.  3  —  Time  history  of  mode  amplitudes  A2to  (along  initial  cloud  density  gradient) 
and  Aifi  for  L  *  10  km  where  Am/t  -  Amkw,„k„  and  kox  «  koy  *  2ir/960  m.  The 
horizontal  solid  lines  are  the  steady  state  amplitudes  Ai,iTH  and  A2,o™  from  the 
theory  of  Chaturvedi  and  Ossakow  f  19791. 
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01CY  ATTN  ADWATE  MAJ  BRUCE  BAUER 
OICY  ATTN  HIT 

OiCY  ATTN  DOK  CHIEF  SCIENTIST 

SAt«0/YA 
P.  0.  BOX  92960 
WORLOWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

OICY  ATTN  YAT  CAPT  L.  BLACK  WELDER 

SAMSO/SK 
P.  0.  BOX  92960 
WORLONAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

OICY  ATTN  SKA  (SPACE  COW  SYSTEMS)  M.  CLAVJN 
SAMSO/W 

NORTON  AFB,  CA  92409 
(MINUTEMAN) 

OiCY  ATTN  WN.  LTC  KEfPCOY 
COTWCER 

ROME  AIR  DEVELOPMENT  CENTER,  AF$C 
HANSCOM  AFB,  MA  01731 

OiCY  ATTN  EEP  A.  LOR ENT ZEN 


DEPARTMENT  OF  ENERGY 


DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.  0.  BOX  3400 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

DEPARTCNT  OF  ENERGY 
LIBRARY  ROOM  G-042 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  FOR  A.  LABOWlTZ 

EGSG,  INC. 

LOS  ALAMOS  DIVISION 

P.  0.  BOX  809 

LOS  ALAMOS,  NM  85544 

OtCY  ATTN  OOC  CON  FOR  j.  BREEDLOVE 

UNlVtRSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.  0.  BOX  808 
LIVERMORE,  CA  94550 

OiCY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 
OiCY  ATTN  OOC  CON  FOR  1-389  R.  OTT 
OICY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 
OICY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 


LOS  ALAMOS  SCIENTIFIC  LABORATORY 

P.  0.  BOX  1663 

LOS  ALAMOS,  W  87545 


OICY 

ATTN  DOC 

CON 

FOR 

J. 

WOLCOTT 

oicy 

ATTN  DOC 

CON 

FOR 

R. 

F.  TASCHEK 

COMMANDER 

OICY 

ATTN  DOC 

CON 

FOR 

E. 

JO*S 

FOREIGN  TECHNOLOGY  01 VI SION,  APSC 

OICY 

ATTN  DOC 

CON 

FOR 

J. 

HALlK 

WRIGMT-PATTERSON  AFB,  OH  45433 

OICY 

ATTN  DOC 

CON 

FOR 

R. 

JEFFRIES 

OICY  ATTN  NICO  LIBRARY 

OICY 

ATTN  DOC 

CON 

FOR 

J. 

ZIW 

OICY  ATTN  ETDP  B.  8 ALLARD 

OICY 

ATTN  DOC 

CON 

FOR 

P. 

KEATON 

OICY 

ATTN  DOC 

CON 

FOR 

D. 

WESTERVELT 

COWAfOER 

RO*€  AIR  DEVELOPMENT  CENTER,  AF$C 

SAfCIA  laboratories 

GRIFF ISS  AFB,  NY  IJ441 

p.  0.  BOX 

SIOO 

OiCY  ATTN  DOC  LIBRARY/TSLD 

AL8uque*que,  wi  ims 

OICY  ATTN  UCSE  V.  COY* 

OICY 

ATTN  DOC 

CON 

FOR 

J. 

HART  IN 

OtCY 

ATTN  DOC 

CON 

FOR 

w. 

BROWN 

SAMSO/SZ 

OICY 

ATTN  DOC 

CON 

FOR 

A. 

TMORNWtW. 

POST  OFFICE  BOX  92960 

OICY 

ATTN  OOC 

CON 

FOR 

T. 

WRIGHT 

WORL>AY  POSTAL  CENTER 

OICY 

ATTN  DOC 

CON 

FOR 

D. 

DAH.GAEN 

los  angeles,  ca  90009 

OICY 

ATTN  DOC 

CON 

FOR 

3i*»l 

(SPACE  DEFENSE  SYSTEMS) 

OICY 

ATTN  DOC 

CON 

FOR 

s*»cii  °*CviCT  31  v 

OiCY  ATTN  SZJ 


SANOiA  LABORATORIES 
Ll vERMORE  ,A8OAAT0«y 
P.  0.  BOX  969 
11 vERMORE,  CA  9  **550 

OICY  AT  Tv  DOC  CON  FOR  a.  ^RPhEy 
OICY  ATTN  DOC  COn  FOR  T.  COCA 

OFFICE  OF  MILITARY  APPLICATION 
DEPART^nT  of  ENERGY 
mAShINGTCN,  J.C.  20545 

Oicv  ATTN  DOC  CON  FOR  0.  GAlE 


OThER  GOv€R^€VT 


CENTRAL  INTELLIGENCE  AGENCY 
ATTN  RD/SI#  Rm  5G48,  hQ  BLDG 
WASHINGTON,  D.C.  20505 

01CY  ATTN  0SI/PS1D  RM  5F  19 

DEPARTMENT  of  commerce 
nation*.  Bureau  of  standards 

WASHINGTON,  D.C.  20254 

(ALL  CORRES:  ATTN  ScC  OFFICER  FOR) 
OICV  ATTN  R.  MOORE 


INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECCmjNICATIONS  &  INFO  AOMIN 
BOULDER,  CO  80503 

OICV  ATTN  A.  JEAN  (UNCLASS  ONLY) 
OICT  ATTN  M.  UTLAUT 
OICV  ATTN  D.  CROf-eiE 
OICV  ATTN  L.  BERRY 

NATIONAL  OCEANIC  (  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 

OICV  ATTN  R.  GRUBB 

OICV  ATTN  AERONOMY  LAB  G.  REID 


AEROSPACE  CORPORATION 
P.  0.  BOX  92957 
LOS  ANGELES,  CA  90009 

OICV  ATTN  I.  GARFUT*EL 
OICV  ATTN  T.  SALMI 
OICV  ATTN  V.  JOSEPHSON 
OICV  ATTN  S.  BONER 
OICV  ATTN  N.  STOCKnELL 
OICV  ATTN  D.  OLSEN 
OICV  ATTN  J.  CARTER 
OICV  ATTN  F .  MORSE 
Oicv  ATTN  SMFA  FOR  PWW 

ANALYTICAL  SYSTEM  ENGI6££R|NG  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

OICV  ATTN  RAOIO  SCIENCES 

BEIKELEY  RESEARCH  ASSOCIATES,  INC. 
P.  0.  BOX  983 
BERKELEY,  CA  9*701 

OICV  ATTN  J.  WOfctHAN 


BOEING  COMPANY,  The 
P.  0.  BOX  3707 
SEATTLE,  wA  98124 

Oicv  ATTN  G.  KEISTER 
OICV  ATTN  D.  MURRAT 
01CY  ATTN  G.  rtALw 
OICV  ATTN  J.  KEWEY 

CALIFORNIA  AT  SAN  DIEGO,  JNN  OF 
IPAPS,  0-019 
LA  JOLLA,  CA  92093 

OICV  ATTN  HENRY  G.  BOOKER 

BROWN  ENGINEERING  COMPANY,  INC. 

Cummings  research  park 
HUNTSVILLE,  Al  35807 

01CY  ATTN  ROMEO  A.  DElIBERtS 

Charles  stark  draper  laboratory,  inc. 

555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 

01CY  ATTN  D.  0.  COX 
01CY  ATTN  j.  p.  GILMORE 

COMPUTER  SCIENCES  CORPORATION 
6565  ARLINGTON  3LVD 
falls  ChuRO,  VA  22046 

01CY  AT T\  H.  fiLAMt 
01CY  ATTN  JOHN  SPOOR 
Q1CY  ATTN  c.  NAIL 

COMSAT  laboratories 
l I  nth i CUM  road 
Clarksburg,  md  20734 

OICT  ATTN  G.  HYDE 
CORNELL  university 

department  of  electrical  engineering 
IThACA,  NY  N850 

OICT  ATTN  D.  T.  PARLEY  jR 

electrospace  systems,  inc. 
box  1359 

RICHARDSON,  TX  75080 

Q1CY  ATTN  H.  LOGSTON 

oicy  attn  security  (pa*,l  Phillips) 

ESl  INC. 

495  JAVA  DRIVE 
SuWvalE,  CA  94066 

OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAf€S  MARSHALL 
OICY  ATTN  C.  W.  PRETTIE 

FORD  AEROSPACE  l  OmJNICATIONS  CORP 
3939  FABIAN  MAY 
PAlO  ALTO,  CA  94303 

OICY  ATTN  J.  T.  MATTINGLEY 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALiEY  FORGE  SPACE  CENTER 
GOOOA<J  BlvD  KING  OF  PRUSSIA 
P.  C.  Box  8555 
Philadelphia,  pa  19101 

OICY  ATTN  M.  H.  BORTNER  SPACE  SCI  LAB 

general  electric  COMPANY 
p.  0.  BOX  1122 
SYRACUSE,  NY  13201 

OICY  ATTN  f.  reibert 

GE«RAl  ELECTRIC  COMPANY 
TEIPO-CENTER  FOR  AOVANCEO  STuOIES 
816  STATE  STREET  (P.O.  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OICY  ATTN  DAS I  AC 
OICY  ATTN  DON  CHANDLER 
OICY  ATTN  TOM  BARRETT 
OICY  ATTN  TIM  STEPHANS 
OICV  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICV  ATTN  B.  GAMBILL 
OICV  ATTN  MACK  STANTON 
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ASSOCIATES 
a.  0.  80*  9*3 ?S 
MARINA  DEl  REv,  CA  902vl 

OICY  AT  *\  cORR£$T  *IlMOR£ 

01CV  ATT\  BRYAN  GABBARD 

01CV  ATTN  *iwU4M  B.  rtAIiinT  JR 

OUT  ATT\  ROBERT  r.  lELEVIER 

01CY  ATTN  WILLIAM  J.  KARZAS 

01CV  ATTN  h.  ORY 

OICY  A-TN  MACDONAlD 

OICY  ATT.\  R .  T^RCO 

RA>£>  CORPORATION,  Th£ 

1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

01CY  ATTN  CULLEN  CRAIN 
01CY  ATTN  ED  8E0R0ZIAN 

RIVERSIDE  RESEARCH  INSTITUTE 
SO  WEST  ENO  AVENUE 
New  YOR<,  NY  10023 

OICY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 
p.  0.  BOX  2351 
LA  JOLLA,  CA  92038 

OICY  ATTN  LEWIS  M.  LlNSCN 
OICY  ATTN  DANIEL  A.  HAM.IN 
OICY  ATTN  D.  SACHS 
OICY  ATTN  E.  A.  STRAKER 

oicy  attn  Curtis  a.  smith 

OICY  ATTN  JACK  MCDOUGALL 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
Sudbury,  ma  o 1776 

OICY  ATTN  BARBARA  ADAMS 


:d:V:  -r->_ ;  . 
^•*00  al 5’*'-^*. 
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1 JNS ,  ;  v; . 

30  MISSION 

i  »= 

Pleasanton, 

;a 

>•*5  vj 

C1CV 

A"r\ 

SC 

SRI  INTERNATIONAL 

333  RAvENS*.OOD  Avtxuc 

ME»*.0  PARK. 

,  CA 

9**  J  25 

OICY 

ATTN 

JUNAwJ  NElLSUN 

OICY 

ATTN 

alan  burns 

OICY 

ATTN 

G.  SMITH 

OICY 

ATTN 

L.  L.  COBS 

OICY 

ATTN 

DAVID  a.  JO-NSON 

OICY 

ATTN 

-ALTER  G.  ChcSNuI 

OICY 

ATTN 

Charles  l.  ring 

OICY 

ATTN 

WALTER  jaye 

OICY 

ATTN 

M.  BARON 

OICY 

ATTN 

RAT  L.  lcAOABRANC 

OICY 

ATTN 

G.  CARPENTER 

OICY 

ATTN 

G.  PRICE 

OICY 

ATTN 

J.  PETERSON 

OICY 

ATTN 

R.  HAKE,  JR. 

OICY 

ATTN 

V.  GONZALES 

OICY 

ATTN 

D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEOFORD,  MA  01730 

OICY  ATTN  W.  P.  BOQUIST 

TRW  DEFENSE  i  SPACE  SYS  GROUP 

one  space  Park 

REDONDO  BEACH,  CA  90278 

OICY  ATTN  R.  K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 


V1SIDYNE,  INC. 

19  third  avenue 

NORTH  WEST  INDUSTRIAL  PARK 

Burlington,  ma  01803 

OICY  ATTN  CHARLES  HU*»HREY 
OICY  ATTN  J.  W.  CARPENTER 


SCIENCE  APPLICATIONS,  INC. 
HUNTSVILLE  DIVISION 

2109  w.  Clinton  avenue 

SUITE  700 

HUNTSVILLE,  At  35805 

OICY  ATTN  DALE  H.  01  VIS 
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